Interacting spin waves in the ferromagnetic Kondo lattice model 
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We present an new approach for the ferromagnetic, three-dimensional, translational-symmetric 
Kondo lattice model which allows us to derive both magnon energies and linewidths (lifetimes) 
and to study the properties of the ferromagnetic phase at finite temperatures. Both "anomalous 
softening" and "anomalous damping" are obtained and discussed. 

Our method consists of mapping the Kondo lattice model onto an effective Heisenberg model by 
means of the "modified RKKY interaction" and the "interpolating self-energy approach". The 
Heisenberg model is approximatively solved by applying the Dyson-Maleev transformation and 
using the "spectral density approach" with a broadened magnon spectral density. 

PACS numbers: 75.30.Mb, 75.50.Pp, 75.30.Ds, 75.10.Jm 



I. INTRODUCTION 

The Kondo lattice model 1 describes the interaction be- 
tween two groups of electrons. One group consists of 
itinerant conduction band electrons which can hop to 
different lattice sites. The other group concerns local- 
ized electrons that couple to a magnetic moment of spin 
S localized at a certain lattice site. Both subsystems can 
perform an intra-atomic interaction with each other while 
neglecting interactions between the itinerant electrons or 
between the localized spins. For non-degenerated band 
electrons in real space, the Hamiltonian reads 

J, 



H — ^2 Tijc\ a c. ](y — — ^ {^z a S^c\ IJ Ci < y + S i c 'c-^a- 

ija ia 
Za = 5 a ] — Sal, ^ = S + S a l + S~5 a l 



(la) 
(lb) 



where c^J represents an annihilition (creation) operator 
for an electron of spin projection a at a lattice site R^. J 
is the Hund's coupling constant and Ty are the hopping 
integrals. Since we are investigating the ferromagnetic 
Kondo lattice model, J > 0. 

The second term in Eq. (flaj) describes an Ising-like 
interaction between the ^-components of the localized 
and the itinerant spin. The third term accounts for the 
spin exchange processes between the two subsystems. 
We will treat the three-dimensional, translational- 
symmetric, infinitely-extended Kondo lattice model. 

The Kondo lattice model is believed to character- 
ize the basic physics of a wide variety of solid state 
materials. 

Magnetic semiconductors, e.g., EuO, are a prominent 
class of substances, which draw notable attention due 
to the "red shift" of the optical absorption edge upon 
cooling from T — Tq to T — OK. One can conclude that 
the coupling constant J is positive and of the order of 
some tenth of eV. In contrast, the magnetic ordering of 
the localized spins is explained via special superexchange 
mechanisms. 



Ferromagnetic local moment metals, such as Gd, are 
another application. An RKKY- (Ruderman and Kittel 2 , 
Kasuya 3 , Yosida^) type interaction is supposed to create 
the ferromagnetic order. The magnetism relies on 
localized 4/ electrons that are shielded from the 4/ 
orbitals of adjacent atoms by other completely filled 
orbitals. On the other side, the conductivity properties 
are determined by itinerant 5d or 6s electrons. 
The discovery of the "colossal magnetoresistance" 
(CMR) and its promising technological application 
motivated a considerable research effort that is related 
to the manganese oxides with perovskite structures 
T^ x D x Mn0 3 (T=La, Pr, Nd; D=Sr, Ca, Ba, Pb). A 
prototype is the well-known compound Lai_a;Ca K Mn03 
which can be obtained by replacing a trivalent La 3+ ion 
with the divalent earth-alkali ion Ca 2+ in La 3+ Mn 3+ 03 
leading to a homogeneous valence mixture of the man- 
ganese ions Mn^ x Mn^+. The three 3d - t 2g electrons 
of Mn 4+ are considered as localized forming a spin of 
S = §. They are coupled to the n — (1 — x) itinerant 
3d — e g electrons per Mn site by a ferromagnetic 
coupling J > 0. J is estimated to be much larger than 
the electronic bandwidth since the manganites are bad 
electrical conductors. 

Many fascinating features of the Kondo lattice model 
can be accredited to the complex correlation between 
the magnetic and electronic subsystem. In this regard, 
one challenging issue represents the "anomalous soften- 
ing" of the spin wave dispersion, which has attracted 
comprehensive interest. The spin wave dispersion 
relation of manganites with high Tq resembles one 
of a simple Heisenberg model with nearest-neighbour 
exchange only>£.i£ However, some manganites with lower 
Tc exhibit apparent deviations from this behaviour, that 
are strongly dependent on the band occupational^ 
Despite extensive theoretical work in this field, the 
softening of the dispersion relation near the boundaries 
of the Brillouin zone still lacks a complete explanation. 
Currently, disorder induced softening has been excluded 
for some materials^ On the other hand, the incorpora- 
tion of an antiferromagnetic super exchange interaction 
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between the Mn ions into the Hamiltonian of the Kondo 
lattice model has been proposed to take into account 
the antiferromagnetic tendencies of the parent material 
LaMn0 3 J£ 

In recent years, unusually large magnon damping at 
the Brillouin zone boundaries and low temperatures 
has come to the fore. This is frequently referred to as 
"anomalous damping."— i Evidence has been found in 
neutron scattering experiments with manganites and 
raised questions concerning the nature of anomalous 
damping and its link to anomalous softening. Besides 
the electron-magnon interaction, some authors speculate 
about a magnon-phonon coupling for certain manganites 
as an origin/ while other authors reject it.— Thus, 
it is of crucial importance to develop new spin wave 
theories for the Kondo lattice model and to study 
whether anomalous damping can be traced back to the 
electron-magnon interaction. 

In this work, we concentrate on the magnetic sub- 
system of the Kondo lattice model. The aim is to 
investigate the dependencies of the energy as well as the 
linewidths of the elementary magnetic excitations called 
spin waves or magnons, respectively. A subsection of the 
paper will treat the anomalous features of the magnon 
spectrum mentioned above and include a discussion of 
the influence of temperature. 

We will introduce a new solution for the Kondo lattice 
model which is as well applicable to the Heisenberg 
model. It goes explicitly beyond standard methods like 
the "random phase approximation," by accounting for 
correlations of higher order. Although non-pertubative, 
it is still controlled in the sense that, in principle, it 
results from the moments of an exact high energy ex- 
pansion. We assume quantum spins, so our theory is not 
restricted to the classical limit of large spin values S> 1. 

The paper is structured as follows: First, we will 
demonstrate how the Kondo lattice model is mapped 
onto a Heisenberg model (Sec. IH|) . Both employed 
theories, the "modified RKKY interaction" and the 
"interpolating self-energy approach", have been already 
successfully applied to the Kondo lattice model for 
various other problemsj 11 ' 12 ^ 3 ' 14 ! 15 They will fix the 
electron-spin interaction. In Sec. IIII1 we will focus on 
the spin-spin interaction by introducing a new solution 
for the Heisenberg model that will not only allow us to 
calculate the energy of the magnetic excitations, but 
also their linewidth. Sec. IIVI will proceed with numerical 
results for the Kondo lattice model that provide insights 
into the properties of its ferromagnetic phase and an 
investigation of the dependence on the intra-atomic 
coupling constant J, the conduction band occupation n, 
and the temperature T . 



II. MAPPING ONTO A HEISENBERG MODEL 

The Kondo lattice model provokes a complex many- 
body problem solvable only in a few limiting cases. 
Hence, we try mapping the Hamiltonian of the Kondo 
lattice model flTa|) onto an effective Heisenberg Hamilto- 
nian in which the conduction band electrons mediate the 
indirect exchange interaction between the localized spins. 
The idea is to use the " modified RKKY interaction^^" 
(mRKKY), wherein the Hamiltonian is averaged in the 
electronic subspace 

H * = -|E (toSfi&Ci,)® +Sr<cU-*) (a) ) • (2) 

The arising expectation value {c ia Ci- a )^ does not vanish 
generally since the spin conservation is valid for the total 
system of the localized spins and the itinerant electrons 
while the averaging is done in the electronic subspace 
only. The expectation values can be calculated by us- 
ing the spectral theorem and the corresponding electron 
Green's functions called "restricted Green's functions" 



(4c lCT ) (s) = ~^f dE f-( E ) Im G T {S) (3a) 

G^ (s) = ((^;4)) (s) (3b) 

(cU-*> W = ~^\ dE f - {E) Im (4a) 

Gr.-W = (( Ci _ CT ; C t £r ))C 8 ), ( 4b ) 



where f—(E) is the Fermi function. After introduc- 
ing the free Green's function for non-interacting elec- 
trons G\j(E), the equations of motion of G™'^ (E) and 
q-<t<t,(s) tjg\ can k e so i ve (j 

G^ s \E) = Gf(E)-±J2 G u- (5) 

i 

G ;p^{E) = -^£ G f- (6) 

I 

■ (-^f((^;i)) (B) +Sf((c la ;c) a ))U>) . 

Now we replace the restricted Green's functions on the 
right-hand sides of Eqs. ([5]) and ^ with their full equiv- 
alents 

G°p is) (E)^G ijr ,(E) = {{c i( ,-,cl)) (7) 
G^ s \E)^G-f°(E)=0. (8) 

G^ (T ' 7 {E) vanishes because of spin conservation. G i j cr (E) 
labels the Green's function of interacting electrons. 
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These solutions are inserted into the averaged Hamilto- 
nian $2§. Eventually, our approach leads to an effective 
Heisenberg Hamiltonian^ 



mRKKY 



Kondo 



(9) 



The effective exchange integrals are functionals of the 
electronic self-energy Y>ij a (E) via the electron Green's 
function Gij a (E) 

Jij = ^ffl dE f - {E) Im Y, G ?h E ) G ^ E )- ( 10 ) 

The replacements (JS|) and © comprise many-body cor- 
relations of higher order than the conventional RKKY 
theory that would be obtained by replacing 



G^ {S) (E)^G^(E), 



-,(0), 



o. 



Hi) 



We use the interpolating self-energy approach 17 (ISA) 
in order to set the electronic self-energy Y, i j <J (E). It is 
derived for the limiting cases of the ferromagnetically- 
ordered semiconductor, the atomic limit, and second or- 
der pertubation theory assuming vanishing band occupa- 
tion n = J2cr( c lcr c icr) ~* 0- An interpolation between the 
limiting cases performed by fitting leading terms in its 
rigorous high-energy expansion provides the result 

Z ijl7 {E)=-^z a (S z )5 ij + (12a) 

| ,P a a Gf(E-\Jz a (S*)) 5 
4 l-6 CT Gf (E-$Jz a {S')) lJ 

a a = 5(5 + 1) - z a (S z )(z a (S z ) + 1), b a = J -. (12b) 

Although derived in the low concentration limit, we 
apply the self-energy (|12a|) to the case of finite band 
occupations n > 0. 

In summary, the problem is reduced to the solu- 
tion of an effective Heisenberg model with exchange 
integrals Jij that depend on the coupling constant 
J, the band occupation n, and the temperature T: 

Jij Jij ( J) ^5 Tj . 



III. THE EFFECTIVE HEISENBERG MODEL 
A. Solution 

It is convenient to transform the spin operators Sj 
of the effective Heisenberg Hamiltonian @ into bosonic 

magnon operators |cii, aj, n.i — aja.ij by means of the 

Dyson-Maleev transformatio n 18 ' 19 ! 20 : 



S+ = V2Sa t , Sr = V2Sa\ (l - ||) 



Sf 



S — m. 



(13a) 
(13b) 



After a Fourier transformation, the Heisenberg Hamilto- 
nian then reads in momentum space 



H 



Hciscnbc 



q' 



(14) 



at4 y , (^q4 ^qi-q3) a qi a q 2 a q3 a q4^qi+q2,q3+q4 



N 4 

qi...q 4 

where fkOq = 25 ( Jo — J q ) stands for the bare energy of a 
free magnon with momentum q and N for the number of 
lattice sites. The second summand of Hum in Eq. (fl4|) 
describes the magnon-magnon interaction and causes the 
existence of finite linewidths and the renormalization of 
the magnon energy. The Dyson-Maleev transformation 
makes it possible to take the complete interaction 
between the magnons into account without making 
approximations that are necessary for other theories, 
e.g., the Holstein-Primakoff transformational. 
At this stage, we need to mention that the transfor- 
mation (|13a|) and (|13b[) leads to unphysical states for 
temperatures near the transition temperature Tq since 
we transform from a Hilbert space which is (25 + 1) 
dimensional into one with infinite dimensions. Neverthe- 
less, according to Dyson, the contributions to the free 
energy from these unphysical states are smaller than 
exp (— ct^S?) where a is a coefficient of order unityji 9 - 
Additionally, Sf and S~ are not Hermitian conjugated 
in the Dyson-Maleev formalism. However, Bar'yakhtar 
et ali2£ showed that the contributions to spin corre- 
lation functions from unphysical states arising from 

the non-Hermiticity are of the order exp ^— where 
k B T* = 5(25 + 1)J . 

The bosonic Heisenberg model (IT41 is solved by 
applying the "spectral density approach." The spectral 
moments of the spectral density Sq(E) are defined by 



(15) 



but they can also be evaluated exactly and independently 
from Eq. (fT5|) by the following relation 



M {n) = 



(16) 



p-fold commutator (n— p)-fold emmutator 

The approach requires a spectral density which is usually 
guessed, e.g., from experiments or theoretical considera- 
tions. Parameters of Sq(E) can be evaluated by a suffi- 
cient set of equations that is derived from the equivalence 
of Eqs. JT5J) and (HBJ). 

In our case, 5 q (i?) represents the magnon spectral den- 
sity which is associated with the average magnon occu- 
pation number by the spectral theorem 



(a q a q ) = (n q ) = ± j dE f+(E)S^E), 



(17) 
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where f+(E) is the Bose function. Since we are inter- 
ested in lifetime effects, we need to fit the first three 

spectral moments M q ™- 2 ^ and use a spectral density of 
finite width. The renormalized magnon energies M7 q and 
their spectral linewidths T q or lifetimes, respectively, 



(18) 



derived by decoupling them using a mean field approxi- 
mation with respect to momentum conservation 



( a q3+q4 -q a q3 a q 2 a qi+q 2 -q) ^qi ,q4 
MF 



(24) 



i a q 3 +q 4 -q a q3 ) ( a q 2 a qi +q2 -q) <^qi ,q4^q4 ,q+ 
^ a cj3+q4 -q a qi +q2 -q) ( a q 2 a q3 ) ^qi ,q4 <^q 2 ,q3 



work as parameters that need to be calculated from the 
set of equations. 

For simplicity and without loss of generality, we want 
to restrict our derivation to a symmetric spectral density 
Sq(tiSlq + E) = Sqihtlq-E)^ This is in agreement with 
neutron scattering experiments^ and other theories^* 
where Sq(E) has the approximate shape of a Lorentzian 



oLorcntzian/ rp\ q 



or a Gaussian, respectively, 



S. 



Gaussian 



(E) 



(19) 



(20) 



One must keep in mind that the Lorentzian must be re- 
stricted to a finite energy interval ensuring a finite second 
spectral moment M q . 
The zeroth spectral moment 

<) = MJ = 1 (21) 

expects a normalized spectral density Sq(E) according 
to Eq. (fTSf . For the first spectral moment we get: 



M q 1 ) = ([[a q ,F D M]_,a q ]_) 

= ^ q M q °) + ^ ^ (Jq + Jq , __ j __ Jq/ _ q)(r 



(22) 



q' 



The result for the second spectral moment is 
[<2 q , Hdm}_ , [Hbm, a q ] _ 



(23) 



1 

W 



y . (Jq3 + Jq_— q Jq4— q)' 

qi...q 4 

(Jq + Jqi+q2— q — Jqi-q — Jq_ — q)' 



2 (a 



q3+q4 -q a q3 a q 2 a qi +12 -q/ °qi ,q 4 



-n)<5n 



,+ 



^ ( a q3+q4 -q a q 2 a< i3 a qi +q2 -q) ^qi ,q4 



(a qi a) a q3 a q4 ) (5 qi + q2 , q3 + c 



A simple ansatz for the unknown, higher expectation val- 
ues in Eq. ([___]), such as (a q3+q4 _ q a q3 a q2 a qi+q2 _ q ), is 



Therewith, the solution is formally completed: 
2 1 



Mj. 



AT2 



(25) 



))("q')("q")- 



^2(J q / + J q "- q '+ q Jq'-q Jq"-q 

q' q" 

(Jq ~^ Jq" Jq'— q Jq"— q') 
(Jq' "T" Jq" Jq'— q Jq"— q)" 

(Jq + Jq'+q"-q — Jq'— q — Jq"-q 



For numerical reasons, we still need to simplify Eq. (J25J) 
to prevent eight-dimensional integrals in expressions like 

pE ql E q « Jq"-q'+ q Jq"(»v)(»V')- This is done by ex- 
ploiting the translational symmetry: 

JO - Jq = Yj Z ^ i 1 ~ 1«) ( 26a ) 
shells i 

7 (0 = I J2 e iqR . (26b) 



R of 
shell i 



A shell is defined by all lattice sites R at the same dis- 
tance |R — R| to an offset lattice site R4£ Thereby 2„ 
denotes the number of all lattice sites of shell i and Ji 
the exchange integral of shell i. The shells are numbered 
and sorted by the size of their radii, i.e. i — 1 stands 
for the nearest neighbours, i — 2 for the next nearest 
neighbours etc. Within the shell concept, the problem- 
atic terms easily factorize^ 



^ J2 7q"-q'-q 1 7 q ''- q2 (^q')(^q") 



N 2 



(27) 



q'q" 



7 ( 
'q 



q' shells m q" 



where we use a notation similar to Dvey-Aharon and 
Fibich 2 ^: 



ZiZj 



(28) 



Nijm gives the number of shell vectors R m that can be 
constructed out of the sum of shell vectors of the shells i 
and j 



R 2 : -|- Ry — R n 



(29) 



By applying the shell concept not only to the expressions 

( 2) _____ 

of the second spectral moment M q in Eq. ([23)) , but also 
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to the first spectral moment Mq 1 ' in Eq. (|22]) . we finally 
obtain after some algebra 



ftOq =2^> J t (l - 7 «) (5 - A,) 

shells z 
shells i,j 



(30a) 
(30b) 



•(l-7 q l) ) (^^(^-^(l-T^))- 
The influence of the shape of Sq(E) on the second spec- 

(2) 

tral moment A/ q is contained in the dimensionless quan- 
tity m q : 



i y x 2 r q S' q (a;r q + M7 q ). 



(31) 



When using a Lorentzian or a Gaussian spectral density, 
rhq is advantageously independent on the momentum q. 
Fij is merely determined by the Af. 

Fij =A -A t - N ij0 A - N Hm (A - A m ). (32) 

shells m>0 

Therewith, the A; remain the only unknown quantities 
in Eq. (|30ap and (|30b[) at a given temperature. Because 
of the definition 



A = 



lEq'W) if i = 

iEq^l-T?) K') if i > 



(33) 



the Aj depend on the spectral density Sq(E) and via the 
Eqs. (|30ap and (|30b|) on each other 



A = A({^}). 



(34) 



Therefore, they have to be self-consistently calculated. 
With the solution satisfying Eq. (|3"4"|) . one can directly 
compute hflq and T q for a given momentum q. 



B. Comparison to experimental data and other 
theories 

In terms of checking our theory for a real system, 
we consider the Heisenberg ferromagnet EuO, whose 
exchange integrals J\ and J2 are known^ According 
to Fig. [1] a good agreement concerning the magnon 
properties M7 q and T q , and the magnetization (S z ) 
can be found between the numerical results of our the- 
ory and the experimental data for a wide range of low 
and intermediate temperatures and for not-too-small 
momenta! 23 ' 26 ! 27 ' 28 At temperatures near Tq, the un- 
physical states cause a wrong hrst order phase transition 
that contradicts the experimental data. 
Results similar to Eq. (|30b[) have been obtained by other 
theories of the Heisenberg model ) 24 ' 29 ' 30 ' 31 but one notes 
differences for a small range of small momenta q — > 0. 
There, the linewidths of our formula (|30b[) turn out to be 
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FIG. 1: (Color online) Temperature dependence of the 
magnon energy hQ q (T) and the linewidth V q (T) for q = |q| = 
0.8A , and of the magnetization (S Z )(T). The calculations 
(lines) are carried out for the parameters of EuO (ref. [26H in 
comparison to experimental data from the Refs. [23j, 12a . and 
1271 (crosses). As in the experimental analysis, a Lorentzian 
(|19[) is used for the magnon spectral density Sq(E). 



too large: 



r 2 q 2 

+0 



other authors 29 propose at least a 



dependence r q q ^ u q 2 



q - H . This discrepancy must be clas- 
sified as a consequence of our approximations. Further- 
more, our results give Tx(T) ~ T 1A while the authors of 
Refs. [1^ and l30l suggest a stronger dependence r q ~ T 3 . 



IV. THE KONDO LATTICE MODEL 

In order to circumvent the problem of too many pos- 
sible parameter combinations to discuss, we have chosen 
three exemplary regions. Small band occupations should 
be a suitable criterion for ferromagnetic semiconductors 
and J = W for manganites, where W is the electronic 
bandwidth. Intermediate J and intermediate n define a 
parameter range with obvious anomalous magnon soft- 
ening and damping. The setting of the main parameters 
is listed in table |TJ 

Equations (|3"3"|) . (|30ap . and (|30b[) . respectively, predict 
that the A vanish and consequently r q — > for T — > OK 
when no magnon- magnon interaction is present. This 
contradicts the results of other theories of the Kondo 
lattice model which give finite linewidths at T = OK 
due to direct contributions to r q by electron-magnon 



interactions. 32 ' 33 



lattice structure 


simple cubic (sc) 


spin value S 


3 
2 


magnon spectral density Sq(E) 


Gaussian (1201) 


conduction band 


s band, bandwidth W = leV, 




tight-binding 



TABLE I: Setting of the main parameters used in the calcu- 
lations for Sec. IIVI 
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0.4 0.6 
J(eV) 



- J=0.1eV n=0.01 

- J=0.7eV 




-1.5 



-0.5 





E(eV) 



0.5 



1.5 



FIG. 2: (Color online) Small band occupation n = 0.01. Ex- 
change integrals Jij of the first five shells at T = 0.1K as a 
function of the coupling constant J. 



FIG. 4: (Color online) Small band occupation n = 0.01. Spin- 
resolved, local electron density of states pf(E) for T — 5K 
and different values of the coupling constant J. 



A. Small band occupation 

This limiting case is implemented in our calculations 
by setting n = 0.01. For all values of J > 0, the exchange 
integrals are positive making ferromagnetism possible 
(Fig. [2]). The growth of the exchange integrals Jij for 
increasing J is accompanied by a corresponding growth 
of both the magnon energies M2 q and the linewidths T q 
(Fig. [3]). For small J, we find that Ml q and T q are 
nearly independent on the momentum q since all ex- 
change integrals Jij are of the same order of magnitude. 
That is why higher exchange integrals J n >i have great 
influence on M7 q and T q . However, for J ss W, the 
magnon energy and linewidth are mainly governed by 
the nearest- neighbour coupling J%, and higher exchange 
integrals can be neglected. Accordingly, hQq and T q 
converge to the common curve shape of a ferromagnetic 
nearest-neighbours Heisenberg model. 
The spin-resolved electron density of states (DOS) pf{E) 
(Fig. [4} features the typical properties of the ISA for the 
case of low temperatures^ For weak couplings J s» 0, 



there is just a small exchange splitting between pf(E) 
and pf(E). For strong couplings, the pf(E) band splits 
into two sub-bands. One is shifted by about +4(5 + 1) 
to larger energies and is built up by electrons that sta- 
bilize their state by permanently absorbing and emit- 
ting magnons (" magnetic polaron" ) . The second band 
at smaller energies is the scattering band for electrons 
that have flipped their spin by emitting a magnon. At 
non-zero temperatures, a high-energy sub-band for f- 
electrons, too, emerges, mainly provoked by thermally- 
excited magnons. 



B. "Strong coupling" regime 

Although the strong coupling regime is often identified 
with the condition J 3> W, we will use this term for the 
situation J = W as well, since it marks a threshold in J 
whereupon no qualitative deviations appear any more in 
the quantities which we regard here. 
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=0.01 


- J= 


=1.0eV 


- J= 
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- J= 


=0.5eV 
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=0.3eV 


- J= 


=0.1 eV 



X M R 

momentum q 




X M R 

momentum q 



FIG. 3: (Color online) Small band occupation n = 0.01. 
Magnon dispersion relation ftfi q and momentum dependence 
of the linewidths T q for T = 0.1K and different values of the 
coupling constant J. 



FIG. 5: (Color online) "Strong coupling" regime J = 1.0eV\ 
Magnon dispersion relation ftfiq and momentum dependence 
of the linewidths F q for T = IK and different values of the 
band occupation n. 
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FIG. 6: (Color online) "Strong coupling" regime J = l.OeV . 
Exchange integrals J,j of the first five shells for T = IK as a 
function of the band occupation n. The broken lines show the 
behaviour according to Eq. (|12a[) with the setting (S z ) = S 
for values of n where ferromagnetism is impossible at T = IK. 
The inset shows the temperature dependence of the exchange 
integrals Jij(T) of the first five shells for n — 0.2 and temper- 
atures up to Tc 



FIG. 7: (Color online) Band occupation n = 0.25. Magnon 
dispersion relation Ml q and momentum dependence of the 
linewidths F q for T = IK and different values of the coupling 
constant J. For couplings below J = 0.45eV, the equation for 
F q (Eq. (J30bj) ) yields no real solution at T — IK, indicating 
strong antiferromagnetic tendencies in the exchange integrals 

Jij ■ 



Starting at small band occupations n and with increasing 
n, the magnon energies hflq grow and the linewidths T q 
decline (mainly at the X point, Fig. This is made 
clear by the fact that the exchange integrals Jij first grow 
because more indirect exchange between the localized 
spins is possible when there are more conduction band 
electrons present (Fig. [6]) . M7 q and T q reach a maximum 
(minimum) at about quarter band filling n ~ 0.5 and 
take the usual curve shape of a ferromagnetic nearest- 
neighbours Heisenberg model. When reaching even larger 
electron densities n, this behaviour is reversed: The ener- 
gies decrease and the linewidths increase (mainly at the 
R point), relying on negative higher exchange integrals 
Jn>i < (antiferromagnetic coupling) and the declin- 
ing nearest-neighbour coupling J\ . For band occupations 
above a critical value of n ~ 0.8, this trend leads to neg- 
ative magnon energies that destabilize the ferromagnetic 
order and prefer antiferromagnetism. The reason why the 
antiferromagnetic state is more favourable for the case of 
half band filling n = 1 can be understood with Pauli's ex- 
clusion principle: The electrons can reduce their energy 
when they virtually hop to adjacent lattice sites, which 
is only possible if there is no electron with the same spin 
present. 

A study of the temperature dependence (inset of Fig. 
[6]) reveals that primarily the nearest-neighbour coupling 
J i is enhanced for rising temperatures, while higher ex- 
change integrals remain temperature independent. 



C. Anomalous Magnon Softening and Damping 

In the case of intermediate couplings J and inter- 
mediate band occupations n, the magnon energies and 
linewidths sensitively depend on both J and n. We in- 



vestigate the situation for different values of J in the 
vicinity of a band filling of n — 0.25 and J = 0.4eV. 
Higher exchange integrals J n >i are comparatively large 
and often negative for J < 0.5eV (Fig. [8]), giving rise to 
distinct deformations of the magnon dispersion relation 
hQq and of the curve shape of T q mainly at the Brillouin 
zone boundaries (Fig. [7]). The strongest modifications of 
the magnon dispersion relation occur around the X point 
along with smaller ones at the R point. Below a critical 
J = 0.28eV, parts of the Brillouin zone evolve where the 
magnon energy becomes negative and for this reason fer- 
romagnetism unstable. The linewidths exhibit deviations 
from the common behaviour of a ferromagnetic nearest- 




0.4 0.6 0.8 

J(eV) 



FIG. 8: (Color online) Band occupation n = 0.25. Exchange 
integrals Jij of the first five shells for T = IK as a function of 
the coupling constant J. The broken lines show the behaviour 
according to Eq. (I12a[) with the setting (S z ) = S for values 
of J where ferromagnetism is impossible at T = IK. 
The inset shows the dependence on the distance |Rj — R, | of 
the exchange integrals Jy for T = IK and different values of 
J. The lines are a guide for the eyes. 
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neighbours Heisenberg model between T and X point and 
at the M point. Compared to X and M point, they are 
unusually small at the R point, which leads to unusually 
long magnon lifetimes. When approaching the critical J, 
the linewidths dramatically increase as expected from a 
Heisenberg model near the transition from the ferromag- 
netic to the paramagnetic state. 

Furthermore, we find distinct long-range oscillations of 
the exchange integrals between ferromagnetic Jy > 
and antifcrromagnetic coupling < qualitatively si- 
miliar to the conventional RKKY theory (inset of Fig. 
0). 

It should be pointed out that all the mentioned effects 
are consequences of solely electron-magnon and magnon- 
magnon interactions. 

When we increase the temperature (Fig. 0, M7 q and 
T q reveal unexpected characteristics. The deviations at 
the X point in relation to the usual Heisenberg model 
are reduced, and the magnon energies grow with increas- 
ing T - even at temperatures near Tq. This behaviour 
radically differs from the usual behaviour observed for a 
Heisenberg model (e.g. Fig. [T|). It relies on the growth 
of the corresponding exchange integrals Jy(T), mainly of 
Ji and J4, favouring the ferromagnetic order. Further- 
more and in contrast to the Sees. IIV Al and IIV Bl and to 
the results for EuO in Sec. IhTbI we observe a relatively 
strong temperature dependence of the linewidths, which 
behave like T X {T) ~ T 2 for J = OAeV (Fig. [SJ. 
Although the magnon density of states (inset of Fig. 
contains the characteristic tight-binding curve shape ow- 
ing to a dominating nearest-neighbour exchange Ji, it 
exhibits deviations at energies E w ZOmeV and low tem- 
peratures due to the deformations in M7 q and T q . The 
rise in temperature firstly leads to larger spectral weight 
at E 30meV and finally washes out the structure be- 
cause of the larger linewidths near To 
Anomalous magnon softening and damping can be 
detected in neutron scattering experiments with 
manganitesi&£&&2i From the theoretical point of view, 
anomalous softening at the X point has been reported 
by other author s 11 ' 35 confirming the parameter range of 
intermediate J and n. A theory that features anoma- 
lous magnon damping has been proposed by Pandey et 
a \3£ Therein, the spin operators Si in the Hamiltonian 
of the Kondo lattice model (fTaj) are fermionized by local- 
ized electrons in atom orbitals. An "inverse-degeneracy 
expansion approach" is applied, describing the diagram- 
matic contributions in powers of the inverse number of or- 
bitals incorporated in the calculations. The authors have 
mentioned clear differences from the common behaviour 
of a Heisenberg model for HQq at X and M point and 
for T q between M and R point, even for large J = QW 
and J — > 00. In the double exchange limit J — > 00, other 
authors 33 have found anomalous softening and damping, 
too. 



30 r 



25- 



> 
B 

J 20 
% 



V A 








- T=50K 


M 


- T=100K 




- T=200K 


JI w.. 





n=0.25 
J=0.4eV 




300 



T(K) 



FIG. 9: (Color online) Band occupation n = 0.25. Tem- 
perature dependence of the magnon energy hQx (T) and the 
linewidths Vx(T) at the X point for J = OAeV. For temper- 
atures below T = 27 'K, the equation for F q (Eq. ()30b|l ) yields 
no real solution, indicating strong antiferromagnetic tenden- 
cies in the exchange integrals Jij. 

The inset shows the magnon density of states p magnon (_E) for 
J = OAeV at different temperatures up to Tc- 



V. CONCLUSIONS 

We have presented an approach for calculating the 
magnon energies M7 q and linewidths r q for the Kondo 
lattice model and examined their dependencies on the 
band occupation n, the coupling constant J, and the 
temperature T. Likewise, our ansatz allows us to study 
other interesting quantities such as the electron and 
magnon density of states p(E) or the exchange integrals 
Jij. We have studied it for small band occupation, the 
case of J = W, and for intermediate J and n where 
the magnon spectrum shows anomalies at the Brillouin 
zone boundaries. These deviations can be explained by 
partial antifcrromagnetic indirect exchange between the 
localized spins as a consequence of electron-magnon and 
magnon-magnon interaction. We have demonstrated 
that the deformations of the magnon dispersion relation 
due to anomalous softening become smaller as the 
temperature rises. As mentioned above, these anomalies 
are caused by electron-magnon and magnon-magnon 
interactions only. This finding could permit a better 
understanding of the origin of similar anomalies in real 
materials. 

Note that our method is also applicable directly to a 
pure Heisenberg model with given exchange integrals Jij . 

When comparing our numerical results with ex- 
perimental data of Lao.7Cao.3Mn03)£i2i2i we state 
differences in hflq and T q . Although it is difficult to 
compare theory and experiment without knowledge of 
the electronic band structure, it can be argued that 
the differences can be ascribed to the Hamiltonian (fla|) 
which we have used to derive our results (|30a[) and 
(|30bp . Namely, it does not involve electron-electron, 
spin-spin, or electron-phonon interactions, even though 
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they are regarded as essential for the manganites. The Acknowledgments 

linewidths calculated by our method turn out to be 

too small, 7,33 which can be explained by the absent 

electron-phonon interaction. Moreover, it has been 

shown that the incorporation of an on-site interaction 

U between the itinerant electrons changes the magnon 

dispersion relation drastically— In order to take these This work was supported by the SFB 668 of the 
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